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7 SUMMARY 

5. 

0 

Reentry trajectories, including computations of convective and 
radiative stagnation-point heat transfer, have been calculated by using 
equations for a point-mass reentry vehicle entering the atmosphere of a 
rotating, oblate earth. Velocity was varied from 26,000 to 45,000 feet 
per second; reentry angle, from the skip limit to -20°; ballistic drag 
parameter, from 50 to 200. Initial altitude was 400,000 feet. Explicit 
results are presented in charts which were computed for an initial lati- 
tude of 38 N and an azimuth of 90° from north. A method is presented 
whereby these results may be made valid for a range of initial latitude 
and azimuth angles. 


INTRODUCTION 


The reentry phase of spacecraft flight is one of the most diffi- 
cult design areas the engineer must consider. Since the heating asso- 
ciated with this phase reaches an appreciably great magnitude, it will, 
in most cases, dictate the design of the reentry configuration. The 
aerodynamic loads or decelerations are also important parameters that 
enter the design considerations of a supercircular velocity reentry 
vehicle. J 


For many space missions it is desirable and possible to use a 

reentry mode for return of the spacecraft to earth; this is 
particularly true for unmanned vehicles. Approximate solutions are 
available for descriptions of these parameters in the ballistic reentry, 
^ee refs. 1, 2, and 3, for example.) Exact solutions of the heating 
equations and equations of motion, however, are more desirable for even 
general analysis of the reentry problem. Such exact solutions are 
presented herein. 
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Trajectories have been computed on an IBM 7090 electronic data 
processing system for a two-body system consisting of a point-mass 
reentry vehicle and a rotating oblate earth. Although the results are 
limited to particular combinations of azimuths and latitudes, they 
should be useful for the Wallops Island and Cape Canaveral ranges in 
that for these latitudes reasonable azimuths may be obtained. 

The velocity range considered was from c 6,000 to ^5,000 feet per 
second. The reentry angle was varied from near the skip limit to -2C P. 
The ballistic drag parameters considered were 50* 100, and 200. The 
results are presented in chart form. 


SYMBOLS 


L 

1 

7 

5 

0 


A 

A,B,C 

C D 

F 

g 

h 

L 

m 

Q 

<1 

R 

*e 

Rp 

r 

t 


reference area, sq ft 

constants 

drag coefficient 

force, lb 
gravity, ft/sec^ 

altitude above reference sphere, ft 
latitude, deg 
mass, slugs 

total heat input, Btu/ft^ 
heating rate, Btu/ft^/sec 
dynamic pressure, ipV a ^, lb/sq ft 

radius , ft 

earth's equatorial radius, 20 , 925,832 ft 
earth's polar radius, 20,854,892 ft 
radial distance, Rj + h, ft 
time, sec or min 



iH t*- ir\ o 
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V velocity, ft/sec 

V s satellite velocity, ft/sec 

W weight, mg, lb 

a azimuth, deg 

y flight-path angle, deg 

A longitude, deg 

p Q sea-level density, 2.3769 x slug/cu ft 

p ambient density, slug/cu ft 

a range angle, deg 

o> e earth's rotational speed, 0. 729211508 x lO -1 *' 

Subscripts: 

a relative to a point on the earth 

0 initial conditions 

c convective 

D geodetic 

e equatorial 

E easterly component 

1 inertial 

N northerly component 

n model nose 

P polar 

G geocentric 

R radiation 

r radial component 


radians/sec 
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l local 

L latitudinal component 

A longitudinal component 

A dot over a symbol denotes differentiation with respect to time. 


METHOD 


The trajectory parameters presented herein were computed on the 
IBM 7090 digital computer utilizing spherical polar differential equa- 
tions of motion for a particle near an oblate rotating earth. The 
computer program was obtained from unpublished data prepared by John N. 
Shoosmith of the NASA Space Task Group. The following equations 
describe the accelerations of the vehicle in the longitudinal, lati- 
tudinal, and radial directions: 


A = 


r cos L 


- 2i , (a + o^)cos L + 2rL(A + c%)sin L 


( 1 ) 


L =1| 


r T 9 m \- 

+ g L - 2rL - r(A + o%) sin L cos L 


( 2 ) 


.. F r . 0 

r = ^ + g r + rL^ + ri 


(A + he) 


x2 


p 

cos L 


( 3 ) 


A more detailed discussion of these equations may be found in refer- 
ence 4 . These equations are not presented exactly as those of refer- 
ence 4 . Since the equations of reference 4 were derived for a sta- 
tionary earth with a rotating atmosphere and equations (l) to (4) 
include a rotating earth, the \ of reference 4 becomes A + (% in 
the rotating-earth situation. The terms and symbols of the equations 
have been slightly, rearranged with no effect on the structure of the 
derivation. The A equation was not completely differentiated in the 
reference report but can be readily shown to yield equation (l). 

The altitude h in terms of the radial distance r and the local 
radius is defined as 


w 


h = r - E Z 
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where 


and 



(5) 


The following equations define the drag forces acting on 
In the three aforementioned directions: 


the particle 



(6) 


Fl -4 c da|j oosfLp - Ifl) + sin(L[) - 1^)1 (7) 

F r = -qCrA^- sin(Lj) - Lq) + cos^ - 1^)1 (6) 


The easterly, northerly, and radial velocity components of the 
or particle are described as 


vehicle 


Va £ “ cos L (9) 
V &N = rL cos (L jj - Lq) - r sin^ - 1^ ( 10 ) 
V a r = sin ( L D - 1^) + r cos^ - Lg) (u) 


The total velocity, relative to the earth, 
flight path is defined as 


of the particle along the 


V 


a ~ 



(12) 


The range over the curved earth may be computed from the following 
equation: & 
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where 


cos 


The flight -path 


a = sin L c 
angle 7 a 


Range = cri 


Ro + *V 


sin L + cos L q cos L 


may be obtained from 


cos(A 0 - A) 

the relation 


(13) 


(14) 


sin ?a ■ Y 1 - (15) J 

7 

Reference 5 presents the equations for the components of gravity of an ^ 

oblate earth and the derivation thereof* Reference 6 contains the 
equations describing the atmospheric properties that define the ARDC 
1959 model atmosphere used herein* The approximate convective stagna- 
tion heating equations were obtained from reference 7 and are 


17,600 / p_y / 2 / Va \3 • 15 

\j^n \pJ \Vs/ 


(16) 


Qc ” 


q^-. dt 


(17) 


The approximate equations for radiation heating under an equilibrium 
condition are 


where 


q R 


AR n 



1.78 



(18) 


0 

tl 

PS 

• •v 

0 

(v a < 25,000) 

6.8; B = 12.5 

(25,000 < V a < 30,000) 

0.003; B = 19.5 

( 30,000 < V a < 35,000) 

20.4; B = 12.5 

(35,000 < V a ) 



and 
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Qr = f Or at (19) 

Reference 3 presents a derivation and comparison analysis for equa- 
tions (l8) and (19) • The radiative heating rates presented herein 
have been computed for the stagnation point of a sphere with a radius 
of 1 foot. In order to correlate these data to a sphere having a nose 
radius of other than 1 foot, the desired heating rate may be obtained 
simply by multiplying the rate presented by the nose radius of the 
specific configuration. 

The convective heating rate may be manipulated in much the same 
way. Again the 1-foot nose radius has been used throughout. For a 
specific sphere the convective heating rate presented must be divided 
by the square root of the specific nose radius. 


RESULTS AND DISCUSSION 


Initial Conditions 

The results of the computations are presented in figures 1 to 13 
for the range of reentry conditions investigated. The velocity range 
considered was from 26,000 feet per second to 45,000 feet per second. 

The ballistic drag parameter W/C^A range was from 50 to 200. The 
reentry angle y Q was varied from near the skip limit to 20° below 
the local horizontal. All reentries were initiated at an altitude of 
400,000 feet. Both initial conditions and results are presented for 
earth-relative rather than for inertial conditions. 

The results presented herein were computed for a specific initial 
latitude and azimuth; however, they are not limited to one particular 
latitude and azimuth. If the inertial velocity, earth-relative velocity, 
and reentry angle are held constant, the initial azimuth necessary to 
validate the results presented herein for a chosen latitude, between 
0° and 38° north, may be computed. It can be shown that the inertial 
velocity Vj in terms of the latitude, earth-relative velocity, azimuth, 
and flight -path angle is 


Vq 2 = V a 2 + <% 2 r 2 cos 2 L + 2o%rV a cos L cos 7 a sin oc^ (20) 


If equation (20) is solved for sin 
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sin cLq 


Vi 2 ~ V a 2 ~ co e 2 r 2 oo s 2 L 
2V a o> e r cos L cos 7 a 


( 21 ) 


In order to determine an azimuth for a given latitude, the inertial 
velocity must be held constant for the initial latitude and for the 
chosen iatitude. Applying equation (20) to the initial conditions 
and substituting the results into equation (21) yields 


sin a* = cos L o sin go + cos 2 Lq 

cos L t 


2V & cos L cos 7 


o) e r cos L 
2V a cos 7 


( 22 ) 


fao Ce the results presented herein were computed for a latitude of 
30 north and azimuth of 90° from north, equation (22) becomes 


sin oa = °llg- 801 + 


0.62096cn e r cn e r cos L 


cos L 2V a cos L cos y 2V a cos y 


(23) 


By choosing a desired latitude, flight-path angle, and earth-relative 
e<luatlon ( 2 ^) ma y be solved for an azimuth necessary to 
v ate these results. Figure 1 shows two curves generated from this 
equation or the azimuth and latitude relationships for which the results 

A L Ca ? be Seen from this chart > reentry angle and velocity 
have littie effect on the azimuth change. The two limits presented are 
for the two most extreme cases of all those discussed. The lower limit 
was calculated for an initial velocity of 45,000 feet per second and a 
reentry angle of 0° (which would not enter but Is shown as a limiting 

S S nm? r ! aS the UPper limit Was ca -*- cu lated for an initial velocity of 
2o > 000 feet per second and a reentry angle of -20°. 


Skip Limits 

results from the shallow-angle computations defined an approxi- 
mate skip limit. The skip limit is defined as that reentry angle which, 
for a particular velocity and W/CpA, results ir. a trajectory where the 
reentry vehicle exits from the atmosphere either to escape or to reenter 
farther downrange and may include additional skips before the final 

Crit ! ri0 l f0r , determ lning skip was a trajectory that reached 
an altitude greater than 400,000 feet after reentry was initiated. 

In figure 2 the initial angle y Q at which skip occurs is plotted 
as a function of the reentry velocity for values of W/CpA of 200 and 50. 


P H NlAO 
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The area below these curves is the region where skip will occur. The 
region above the curves Is the region where direct or single pass 
reentry may be obtained. 

Since the accuracy of the skip limit is determined by the number 
of computations made, no attempt has been made to determine precisely 
the skip limit in the present analysis. In determining the skip limit 

no points were computed at intervals less than 1° in y Q . The accuracy 

of the limit is, therefore, considered to be no better than ±-°. There- 

L fore, the data in the charts should be used with caution for reentry 

j angles close to the skip limit defined in figure 2. 

5 

) 

Trajectory Profiles 

Shown in figure 3 are plots of the variation of altitude with range 
for initial velocities from 26,000 to 45,000 feet per second. Profiles 
are given for a range of initial flight-path angles from -20° to a 
shallow angle close to the skip limit. Also shown are time in seconds 
from the time of initiation of reentry and velocity at various points 
along the trajectory. These charts are for a value of w/Cj)A of 100. 

Note that, at the steeper reentry angles, the profiles are essentially 
straight lines down to fairly low altitudes and velocities as predicted 
in reference 2 whereas, for shallow angles, curvature becomes fairly 
pronounced at relatively high altitudes and velocities. Figure 3(d) 

shows that, for V D = 36,000, 7q = -6, and ^ = 100, the reentry 

angle is close to the skip limit j in fact, for a portion of the tra- 
jectory, the reentry vehicle gains altitude. This type of profile is 
typical for trajectories near the skip limit. 

Figure 4 provides a picture of the effect of the variation of 
W/CpA on the altitude-range profile for specific velocities of 26,000, 

36,000, and 45,000 feet per second. The three charts of figure 4 point 
out the marked effect of W/C^ for shallow reentries. As shown in 

figures 4(b) and 4(c) the skip may be determined by w/CjjA. These 
charts also point out the relative insignificance of W/CpA variations 

on the altitude -range profiles at the steeper entry angles. The charts 
of figures 5 and 6 present curves for the variation of altitude with 
velocity for initial velocities of 26,000 to 45,000 feet per second. 
Figure 5 also presents a limited variation of the initial entry angle 
and ballistic parameter W/C-qA for correlation or comparison purposes. 

In figure 6 the midrange velocity of 36,000 feet per second is con- 
sidered. In figure 6(a) the entry angle y Q was varied from -6° 
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to -20°. The drag parameter is varied for the two extremes. Fig- 
ure 6(b) is simply a comparison of the variation of W/CpA. for a velocity 
of 56,000 feet per second and y Q = -10°. As can be noted from these 
charts, the higher W/C-pA flight will dip considerably deeper into the 
velocity altitude corridor than the lower W/CpA flight. 


Range, Time, and Deceleration 

The range from initiation of entry to impact is shown in figure 7 
as a function of entry angle and entry velocity for the three drag 
parameters considered. It can be seen that, a; the skip limit is 
approached, the range curves become asymptotic to the entry angle at 
which skip occurs. Therefore, extreme caution must be exercised in 
using these charts for ranges greater than 1,000 nautical miles. 

A measure of the total time of flight is shown in figure 8. The 
time of flight shall be defined as the time from initiation of reentry 
(that is, 400,000 feet) to impact. The curves for all velocities appear 
to cross near a single point. This apparent crossover is not intended to 
mean that at this point the time of flight is Invariant with velocity, but 
that these curves have been arbitrarily faired through this point. In 
order to define the crossover more accurately a large number of computa- 
tions are necessary. Note that, at shallow reentry angles, the times 
are becoming very large and therefore the trajectory conditions are near 
the skip limit . 

The maximum values of deceleration (along the flight path) are 
shown in figure 9. The curves here again appear to cross at one reentry 
angle. This result can also be attributed to an arbitrary fairing, as 
described above. Reference 2 presents a methol whereby the approximate 
max imu m deceleration may be estimated for steep reentry angles. This 
method provided data that agreed well with that; presented for the 
steeper entry angles (that is, y 0 = -15° to y Q = -20°) but could 
not be compared with more shallow entry angles. 


Heating Rates 

Shown in figure 10 are total convective stagnation-point heat 
inputs (from t = 0 to Impact) plotted as a function of the maximum 
stagnation-point convective heating rate. Convective heating is shown 
for the range of velocities, reentry angle, and W/CpA considered. 

Note that the data are for a spherical nose having a 1-foot radius. 
Data for other radii can be found as described in the section "Method" 
or for other than a sphere by applying appropriate conversions. Fig- 
ure 11 presents the same type of information except that the radiative 
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heating has been added. The total stagnation-point heat inputs are 
plotted as a function of the maximum stagnation-point heating rates. 

Time histories of the stagnation -point heating rates, with integrated 
heating indicated, are shown for a variety of Initial reentry conditions 
in figure 12. 

Shown in figure 13 are the maximum convective and radiative 
stagnation-point heat -transfer rates as a function of initial reentry 
angle for the range of velocity and V/C])A considered. 

As mentioned, the radiative heat-transfer data of figures 11 and 13 
are equilibrium values. The theory for this type of heat transfer is 
fairly well defined. (See ref. 3*) There is some indication (ref. 8, 
for example) that the radiative heat transfer encountered in reentry 
may well be nonequilibrium; however, this type of heat transfer is not 
well understood at present and is often expressed as a multiple of the 
equilibrium value. Therefore, the rates presented can be considered 
as a lower limit of the radiative heating. 


CONCLUDING REMARKS 


Charts have been presented for determination of reentry trajectory 
information, including heating data, for a wide variety of initial 
reentry conditions for velocities of 26,000 to 45,000 feet per second. 
The data are from machine- computed exact solutions of the equations 
of motion. These charts should be useful in preliminary design of 
vehicles entering in a ballistic mode. Although they were computed 
for an initial latitude and azimuth, they may be applied to other 
combinations of latitude and azimuth angles. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Air Force Base, Va. , August 9 y 19^1 • 
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Figure 2.- Skip limits. h Q = 400,000 feet; accuracy = ±0.25°. 
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Figure 3- _ Continued. 
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(f) V Q = 45,000 feet pei second. 
Figure 3 .- Concluded. 
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Figure Altitude-range profile 




Figure 4. - Continued.. 
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Figure 6.- Concluded. 
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Figure J.- Range to impact. 
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Figure J.- Continued 
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Figure 8.- Time to slow to a speed of 10,000 feet per second. 
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Figure 8.- Continued. 
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Figure 8.- Concluded 
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Figure 9- _ Deceleration 
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Figure 9-~ Continued 
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Figure 9* - Concluded. 
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Figure 10.- Total convective heat input as a function of the maximum 
convective stagnation- point heating rate. 
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Figure 10.- Continue! 
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Figure 10.- Concluded. 
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Figure 11.- Total heat input as a function of the maximum total 
stagnation-point heating rate. 
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Figure 11.- Continued 
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Figure 11. - Continued. 
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Figure 11. - Concluded. 


Convective heating rate, q Btu/ftysec 



Figure 12.- Stagnation-point heating rate as a function of time. 
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Figure 12.- Continued. 
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(c) V Q = 56,000 feet per second; y Q = -10°. 


Figure 12.- Continued. 
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Figure 12.- Continued 
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Figure 13*- Maximum heating rates. 
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Figure 13.- Continued. 
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Figure 13.- Concluded. 


